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Abstract —It has been shown that cooperative localization is 
capable of improving both the positioning accuracy and coverage 
in scenarios where the global positioning system (GPS) has 
a poor performance. However, due to its potentially excessive 
computational complexity, at the time of writing the application 
of cooperative localization remains limited in practice. In this 
paper, we address the efficient cooperative positioning problem in 
wireless sensor networks. A space-time hierarchical-graph based 
scheme exhibiting fast convergence is proposed for localizing the 
agent nodes. In contrast to conventional methods, agent nodes 
are divided into different layers with the aid of the space-time 
hierarchical-model and their positions are estimated gradually. 
In particular, an information propagation rule is conceived upon 
considering the quality of positional information. According to 
the rule, the information always propagates from the upper layers 
to a certain lower layer and the message passing process is further 
optimized at each layer. Hence, the potential error propagation 
can be mitigated. Additionally, both position estimation and 
position broadcasting are carried out by the sensor nodes. 
Furthermore, a sensor activation mechanism is conceived, which 
is capable of significantly reducing both the energy consumption 
and the network traffic overhead incurred by the localization 
process. The analytical and numerical results provided demon¬ 
strate the superiority of our space-time hierarchical-graph based 
cooperative localization scheme over the benchmarking schemes 
considered. 

Index Terms —Cooperative localization, space-time 

hierarchical-graph, information propagation, energy 
management, activation mechanism. 

1. Introduction 

Having accurate positional information is a vital require¬ 
ment in various applications, such as traffic control, security 
tracking, search and rescue operations, medical assistance and 
the Internet of Things Q'Pl- This wide range of applications 
has inspired growing interests in localization in wireless sensor 
networks (WSNs) |[8|-|T0|. In a WSN, it is usually assumed 
that anchor nodes whose positions are known can be employed 
for locating the agent nodes. 

In general, the existing positioning techniques may be clas¬ 
sified into two categories: non-cooperative methods and coop- 
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erative methods. For the non-cooperative methods (ID-(E)> 
an agent node is localized by measuring the distances between 
itself and the anchors, whose positions are known. In contrast 
to the non-cooperative methods, cooperation among the agents 
is capable of improving positioning accuracy, wherein an 
agent measures the distances from itself to both the anchors 
and the other agents (HI- Recently, cooperative localization 
has attracted significant research interests |T5|-|[26|. Most of 
the state-of-the-art cooperative localization methods are based 
on the point estimation philosophy, as exemplified by the 
maximum likelihood method the non-linear least squares 
method (H) the two-step linear least squares method fT^ 
and the geometric optimization methods Gz) However, these 
approaches estimate the positions of all the agent nodes in a 
“hard-decision” manner, and they do not take the statistical 
properties of the estimation error into further consideration. 
Therefore, these schemes suffer from the lack of a statistical 
interpretation that could be useful for developing new methods 
exhibiting an improved performance. 

Compared to the above-mentioned point estimation meth¬ 
ods, Bayesian approaches consider the position of an agent 
node as a random variable, and their goal is to estimate the a 
posteriori marginal probability density functions (PDF) of all 
agents’ positions, which is regarded as a localization inference 
problem. Nevertheless, the marginalization of the joint a 
posteriori PDF of all the agent nodes imposes a potentially 
excessive computational complexity, especially in large-scale 
networks. As a remedy, graphical model based methods have 
been proposed in |T5| , p7| for simplifying the localiza¬ 
tion inference problem |T5| , |T^ , pQ| . The message-passing 
methods constitute the class of “best-known” methods which 
make inference on the graph, resulting in an approximately 
optimal solution p5| , p7|. In particular, the belief propagation 
(BP) algorithm |28| represents an efficient message¬ 

passing method of calculating the marginal distribution based 
on probabilistic graphical models, and it can be applied to 
cooperative localization in WSNs |T5| , | [29| . It is noted that the 
computational complexity of the BP algorithm is proportional 
to the total number of links on the graph, thus it is much more 
efficient than the naive “brute-force” method, which imposes 
an exponentially increasing computational complexity when 
computing the marginal probabilities. However, the standard 
BP is unable to resolve a non-Gaussian uncertainty, which 
is often encountered in a practical localization scenario. In 
contrast to the standard BP, the nonparametric belief prop¬ 
agation (NBP) is a sample-based algorithm, which can be 
readily applied in non-Gaussian inference problems p9| . 
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Therefore, the NBP is more suitable for localization. Having 
said this, an impediment of NBP (and also of BP) is that for 
graphs exhibiting loops, there are no guarantees concerning 
the quality of the marginal beliefs as well as the conver¬ 
gence of the algorithm. Moreover, the employment of NBP 
is impractical in dense WSNs, because directly employing 
all the links of the graph for localization without any pre¬ 
selection leads to a high computational complexity, a high 
network traffic, as well as an increased energy consumption. 
In order to mitigate the infiuence of loops in the graph, a 
generalized belief propagation (GBP) algorithm was proposed 
in p0| , p0| . However, GBP also remains unsuitable for large- 
scale WSNs owing to its excessive computational complexity. 
In addition, although the NBP defined over spanning trees 
with breadth first search (NBP-BFS) | [3T| , the tree-reweighted 
belief propagation (TRW-BP) pJl and the NBP defined over 
the minimum spanning tree (NBP-MIN) p2| techniques are 
potentially capable of outperforming NBP when dealing with 
loops, they also have their particular limitations. To elaborate 
a little further, NBP-BFS exhibits poor performance in WSNs 
having sparse nodes, NBP-MIN is less attractive than the 
standard NBP in WSNs having short transmission radius, 
while TRW-BP is unsuitable for a distributed implementation, 
since it requires the knowledge of the PDFs over all spanning 
trees. 

On the other hand, in order to facilitate a more practical 
implementation that imposes both a reduced network traffic 
and a low computational complexity, a technique referred to 
as information censoring has been considered in p2|-p4|, 
which results in carefully controlled information propagation. 
In l l^ , l l^ , the Cramer-Rao bound (CRB) was considered as 
a criterion for removing any unreliable positional information. 
Nevertheless, it is possible that some censored information 
remains helpful to particular agent nodes. Therefore, a more 
judicious control of information propagation (message pass¬ 
ing) is required for improving the positioning quality and 
for reducing the computational complexity. The bootstrap 
percolation approach of p^ , p4| is capable of controlling 
error propagation, hence it has attracted substantial attention 
in recent years. It was shown in p6| that bootstrap percolation 
relying on a hard threshold is capable of carrying out iterative 
localization p5| . However, some agent nodes that do not 
satisfy the above hard threshold cannot be localized with the 
aid of this approach. In our previous work 1^ , we have 
shown that for non-Bayesian cooperative positioning (e.g. 
point estimation), soft threshold aided bootstrap percolation 
is capable of improving the localization performance, despite 
reducing both the computational complexity and the energy 
consumption. 

In this paper, we aim to extend the bootstrap percolation 
approach for solving the Bayesian localization problem. To 
achieve this goal, a novel NBP-based cooperative localization 
scheme is proposed, which is capable of reducing the network 
traffic, whilst accelerating the convergence of the positioning 
algorithm. More specifically, a space-time hierarchical-graph 
based description of the network is proposed, which relies 
on the bootstrap percolation philosophy for facilitating the 
implementation of NBP-based cooperative localization, where 


the agent nodes are divided into different layers, and the 
localization process is carried out in a layer-by-layer manner. 
In particular, based on the framework of the space-time 
hierarchical-graph description of the network, a beneficial 
information propagation rule is conceived for controlling the 
information propagation. This rule includes a joint design 
of inter-layer message passing and intra-layer message pass¬ 
ing. Among different layers, the messages always propagate 
from the upper layers to a certain lower layer. However, 
the messages originating from lower layers, which contain 
the erroneous information, are not allowed to propagate to 
the upper layers. As a result, the impact of any erroneous 
information can be mitigated. Within a single layer, the mes¬ 
sages are allowed to be exchanged among intra-layer sensor 
nodes, when there are insufficient reference nodes (including 
both agent and anchor nodes) in connection. By employing 
the proposed information propagation rule, the positioning 
accuracy can be significantly improved. Additionally, in order 
to manage the energy consumption and network traffic more 
effectively, both position estimation and position broadcasting 
are supported by the sensor nodes. In contrast to traditional 
cooperative methods, we put our emphasis on the behavior 
constraints of sensor nodes by employing a specific activation 
mechanism. Note that the activated agent nodes are capable of 
estimating their positions and of broadcasting their positional 
information, whereas the inactive agent nodes remain silent in 
order to save energy and to reduce the network traffic. Our 
main contributions in this paper are summarized as follows. 

1) An NBP-based cooperative localization scheme, which 
exploits the proposed space-time hierarchical-graph, is 
conceived for WSNs. By formulating the localization 
inference problem subject to additional constraints, a 
reliable solution can be guaranteed. Provided that we 
stipulate the constraint that there should be at least 
three reference nodes in support of position estimation, 
reliable positional information may be gleaned. Addi¬ 
tionally, this constraint ensures that the feasible region 
is shrunken into a small region surrounding the optimal 
solution. 

2) Inspired by the bootstrap percolation strategy, a space- 
time hierarchical-graph is formulated, which is capable 
of mitigating the infiuence of loops contained in the 
graph. The soft connection constraint invoked by our 
strategy ensures that most of the agent nodes share a 
sufficiently high number of reference nodes for their 
position estimation action. Compared to the traditional 
spanning tree method of | [3T| that does not rely on a 
sufficiently high number of reference nodes, the pro¬ 
posed space-time hierarchical-graph is characterized by 
multiple layers emerging from invoking the bootstrap 
percolation strategy. Therefore, the positioning result 
becomes more accurate. 

3) An information propagation rule designed for our space- 
time hierarchical-graph is investigated, which relies on 
the joint design of the inter-layer as well as intra-layer 
message passing. The inter-layer message passing is 
capable of mitigating the error propagation, while the 


ACCEPTED TO APPEAR ON IEEE TRANSACTIONS ON SIGNAL PROCESSING, SEPT 2015 


3 


TABLE I 

List of Major Variables/Notation Used in This Paper 


Variable/N otation 

Definition 

Xz 

Two-dimensional position of agent i, i ^ I •= [1? • • • ? A^] 


Two-dimensional position of anchor j, 
j e J :=[N + 1,...,N + M] 

dyi 

Noisy measurement of the distance between node v and node 

d(x„,Xi) 

Actual distance between node v and node i 

Tiyi 

Noise of the distance measurement between node v and node 

/j 

R 

L 

Transmission radius 

Si 

Set of neighboring nodes of node i 

P(xi) 

Prior distribution of the position of node i 

Q 

A clique which is a subset of V 

V’i(Xi) 

The single-node potential of node i 

V>t.i(x„,Xi) 

The pairwise potential between node i and v 

bU^i) 

Belief (or marginal distribution) on the position of node i at 
time slot t 


Message from node v to node i at time slot t 

7(xi) 

Redundancy at node i to represent the diversity of two 
arbitrary distributions 

Ti 

A set containing the selected sensor nodes for information 
fusion, hence it is a subset of Si. 

Ai 

The set of agents that are activated on layer 1. 

V} 

The set of candidate reference nodes on layer 1 

L 

The number of layers in the hierarchical NBP algorithm 

K 

The number of weighted samples of NBP 


intra-layer message passing ensures that all beneficial 
information can indeed be fully exploited. Additionally, 
with the aid of a beneficial activation mechanism, both 
efficient information control and energy/traffic manage¬ 
ment can be achieved by relying on this rule. 

4) The computational complexity and network traffic of 
the proposed scheme is analyzed. Since both of them 
are related to the number of actively participating links 
of the graph, we adopt the average number of links 
used in the entire network as the metric of compar¬ 
ing the complexity of different NBP algorithms. We 
demonstrate that the proposed scheme is capable of 
significantly reducing the computational complexity and 
network traffic by censoring the unreliable messages 
emerging from the lower layers. 

The rest of this paper is organized as follows. In Section II, 
we describe the system model and formulate the problem 
considered. In Section III, we characterize the localization 
inference problem on a graphical model and propose a novel 
space-time hierarchical-graph based scheme for cooperative 
localization. We present both the computational complexity 
and network traffic analysis of the proposed algorithm in 
Section IV, and characterize the positioning performance in 
Section V. Finally, our conclusions are offered in Section VI. 

Notations'. A/'(/i, cr^) denotes a normal distribution with 
expectation of ji and standard deviation of cr; | • | denotes 
the cardinality of a set; || • || represents the L2-norm of a 
vector; E{-} stands for the expectation operation; 0 denotes 


an empty set; A\B = {x\x G A and x ^ B} denotes the 
relative complement of B in A', H(yi) is the entropy of x; 
I{P} is the indicator function, which has a value of one when 
P is true and zero otherwise; Cov['] denotes the covariance; 
and I is the identity matrix. Table I summarizes some of the 
major variables used in this paper. 

II. System Model 

We consider a two-dimensional WSN consisting of two 
types of sensor nodes: N agents and M anchors. The agents 
have unknown positions at x^, i G / N], and x^ is 

a two-dimensional random variable representing the position 
of agent node i. Similarly, the anchors have known positions 
at xj, j e J := {N + 1,..., A/' -1- M}. Thus, all sensor node 
indices are included in V = I\jJ. We assume that agent node 
i acquires a noisy measurement dyi, which is the estimated 
distance from sensor node v (v can be either an agent or an 
anchor) to agent node i during each time slot. Then, we have 

dyi = d(x^, Xi) -f riyi, G V, i G /, (1) 

where d{^y^^i) = ||x^ ~ ^i\\ is the actual distance between 
node V as well as node i, and Uyi is the measurement noise, 
which may obey log-normal, Gaussian or any other appropriate 
distribution measured in various deployment districts. It is 
worth noting that NBP is applicable to all the possible distri¬ 
butions, since NBP is a sample-based nonparametric method. 
Without loss of generality, we assume that riyi ^ A/’(0,cr^J 
GD and the probability of the event that node v and node i 
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can detect each other is given by 


Pd{Xv,Xi) 


1, ||x„ - Xj|| < R, 

0, otherwise, 


( 2 ) 


where R represents the transmission radius. Note that the 
pair of nodes v and i involved in the detection have been 
connected with each other via an edge, which is denoted 
by satisfying (v^i) G E, where ^ is a set contain¬ 

ing all the edges. Then, for each sensor node i e I, the 
set of neighboring nodes connected to node i is given by 
Si = {i;| ||x^ — x^ll < R, V G V}. Additionally, we assume 
that the position of a sensor node v has an a priori distribution 
of p(x^). For an agent, p(x^) is a uniform distribution 
over the entire area, while that of an anchor is the Dirac 
Delta function at its true position. Furthermore, we assume 
that a) the positions of the agent nodes are independent a 
priori, i.e. we have p(xi,...,x^r) = t)) the 

relative positions during arbitrary time slots are conditionally 

independent and depend only on the two nodes involved, i.e. 

N 

we have p({(i^J|xi,... ,XAr) = [I 11 x^), where 

i=iseSi 

{dvi} denotes the set of measurements concerning the distance 
from sensor node i to sensor node v, i e v e V, i v. 
Then, the joint PDF of the positions of all agents can be 
characterized as 


p(xi,...,XAr,{dm}) 

= n Pd(x^,xi) n p(^it;iix^,xi) n p(x^)- o) 

{v,i)^E {v,i)EE vClV 

Our objective is to compute (or approximate) the a posteriori 
marginal PDF p(yii\{dvi}) for each agent node i. Then, either 
the minimum mean square error (MMSE) estimator or the 
maximum a posteriori (MAP) estimator can be employed 
for completing the location estimation process. In this paper, 
the MMSE estimator that has a relatively low computational 
complexity is invoked, since ideally the position inference 
process should impose a low computational complexity and 
a low additional network traffic, which is also the objective 
of our proposed scheme to be presented in the forthcoming 
sections. 


III. The Proposed Space-Time Hierarchical 
Localization Scheme 

In this section, we first present the statistical framework 
invoked for sensor localization relying on a graph. By formu¬ 
lating the localization inference problem subject to additional 
constraints, finding an accurate solution is guaranteed by the 
proposed scheme. Then, upon exploiting the idea of bootstrap 
percolation subject to soft connection constraints, a space- 
time hierarchical graphical model is formulated. Based on 
this model, we conceive an information propagation rule 
defined over the space-time hierarchical-graph, which requires 
a joint design of the inter-layer message-passing and intra¬ 
layer message-passing. Finally, a new NBP-based localization 
algorithm is conceived, which carries out the localization in¬ 
ference based on the proposed space-time hierarchical model. 



Fig. 1. An example of cooperative localization, where A, B, C and D are 
the agent nodes, while 1, 2, 3, 4 and 5 are the anchor nodes, (a) represents 
an undirected graph, (b) represents a bidirected graph, (c) represents an 
unidirectional graph. 


A. Inference on Graph 

A graphical model is a probabilistic model which employs 
a graph for denoting the conditional independence structure of 
N random variables. Typically, there are two types of graph¬ 
ical models, namely the directed graphical models (Bayesian 
networks) pS] and the undirected graphical models (Markov 
networks) |15| , |T^ , pQ| . As shown in Fig. cooperative 
localization can be defined on either an undirected graphical 
model or a directed graphical model (including unidirectional 
graphs and bidirected graphs). A factor graph (EG) is capable 
of providing a unified interpretation of these graphical models 
and lends itself to a low-complexity implementation of dis¬ 
tributed localization. 

An EG is a bipartite bidirected graph that illustrates how a 
complicated global function relying on many variable factors 
may be simplified into the product of several simple local 
functions, where each local function has fewer variables. More 
specifically, an EG has a variable node representing each 
variable, a function node for each local function, and an edge 
connecting a variable node to a function node if and only if the 
variable is an argument of the local function. The relationship 
between the EG and the joint distribution of all variables 
can be quantified using the potential functions f)Q, which are 
defined on each “clique” Q of the graph. Explicitly, a clique 
is a subset of the complete set V of variable nodes, and every 
two vertices in the clique are connected if the clique contains 
more than one vertices. The joint PDF of all agent positions 
is then formulated as: 

p(xi,...,Xiv) (X jq : q e Q}). (4) 

clique Q 

For the distance-based localization in our system model, the 
joint a posteriori distribution of the sensor positions can be 
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Fig. 2. An FG of the example given in Fig. 1(b). Variable nodes are 
represented by circles, function nodes are represented by squares. Edges that 
connect the agents and anchors are represented by solid lines. Edges that 
connect two agents are represented by dashed lines. Every two variable nodes 
exchange information according to the rule specified by a function node. 


expressed by the potential functions defined over the set V of 
variable nodes as follows 

p(xi,...,XAr|{(l„i}) oc P[t/)„(x„) P[ V)„i(x„,Xi), f e V, 

vGV (v,i)EE 

(5) 

where 'ipv(^v) represents the single-node potential, which is 
defined as the a priori information concerning the position 
of sensor node v, i.e. we have 'ipy{'Ky) = p(x^). Additionally, 
x^) is the pairwise potential between node v and node 


i, which is defined as '0^i(x^,Xi) = pd{'y^vi^i)p{dy 


i)- 


Having defined a statistical frame work for sensor localiza¬ 
tion relying on a graph, we can now estimate the locations of 
agents by applying the classic message-passing algorithms for 
approximating the a posteriori marginals of the agent posi¬ 
tions. We consider the sensor nodes in a WSN, as exemplified 
in Fig. 1(b), as the variable nodes in an FG. Correspondingly, 
the FG function nodes compute the message between a pair of 
sensor nodes in a manner as shown in Fig. 2. The sensor nodes 
estimate their beliefs (or a posteriori marginals) by collecting 
all incoming messages from the neighboring nodes during each 
time slot. Note that the belief of node i during time slot t 
(i.e. in the t-th iteration) is proportional to the product of the 
messages merging into node i and to the local evidence which 
is determined by the single-node potential V^-(xi). Then, the 
updating rule of the belief at time slot t is formulated as 


B|(xj) (X P[ Mliixi), 


( 6 ) 


seSi 


where is the message sent from node s to node i, 

s e Si. The message Ml-((Ki) is computed by the function 
nodes exemplified in Fig. [^according to the message updating 
rule of 


/ 


oeSs\i 


B 


>t-l 






(7) 


The available information is then processed among the 
neighboring nodes according to I© and 0- Each variable 
node modifies its own positional information by fusing the 
local evidence and the received messages. Provided that the 
belief of node i converges, the variable node i becomes capable 
of accurately estimating its position after a few iterations. 
However, convergence may not be reached in those graphs 
which have loops. 

To evaluate the uncertainty of the information propagated 


between the variable nodes and function nodes, we intro¬ 
duce the notion of complementary entropy-ratio (CER) for 
representing the deviation between two arbitrary distributions. 
Explicitly, the CER is defined as 


7(xi) = 1 


JT(x^|{x5, J^J)’ 


5 G Si, 


( 8 ) 


where i7(x*) stands for the entropy of the true distribution of 
the position of node i, while i7(x^|{x5, dg^}) is the entropy 
of the approximate PDF of x^ estimated by the neighboring 
node s. Note that the location information in our approach is 
approximated by a nonparametric representation. As a result, 
it is impossible for the entropy and for the conditional entropy 
to become negative. It is plausible that a lower CER represents 
a more accurate approximate PDF. Before the iterative mes¬ 
sage passing commences, the CER of anchor nodes is zero, 
while that of the agent nodes is one. After several message 
exchanges, a fixed point of convergence can be acquired if 
all messages are consistent and the number of messages is no 
less than three, which implies that the principle of trilateration 
has to be satisfied. Thus, we can formulate our localization 
inference problem subject to additional constraints for the sake 
of achieving a reliable location estimation as follows: 


min 

Ei 

S.t. 


. ^ -1 ^( 4 ) 

7 (xj) = mm 1- - -. 

iJ(xi|{x„,d„i}) 

7(x«) < S, 

Ti C S'*, \Ti\ > 3, 


u e Ti, (9) 


where Ti contains the reference nodes selected for information 
fusion, hence it is a subset of Sp, furthermore, ^ is a small 
value, which indicates that the reference node employed 
provides reliable positional information. Because we have 
7 (x^i) ^0, u ^ J^i, \J^i\ ^3, and Ti C Si, node % 
has three well-positioned neighboring nodes. According to 
the principle of trilateration (B), the “just-sufficiently 
accurate” position of node i can be readily obtained for 
the next-stage process to refine it. Based on the theory of 
NBP ITs) , p(xi|{x^, dui}) (i.e., the approximate PDF of 
node i estimated by the neighboring nodes u e Ti) closely 
resembles p(x*) (i.e., the true distribution of the position of 
node i). 


B. Hierarchical Graphical Model 

In this subsection, a hierarchical graphical model consid¬ 
ering the quality of the positional information of the sensor 
nodes is presented. For notational convenience, we further 
define the set of anchor indices as Tlo = {j\j G J}, and 
treat the set of all the anchor nodes as the root layer in 
the hierarchical graphical model. Additionally, both position 
estimation and position broadcasting are carried out by the 
sensor nodes. In order to control the behaviors of the sensor 
nodes, either the active or the inactive state is assigned to 
each sensor node. Note that only the active nodes are capable 
of estimating the positions and broadcasting messages. In the 
initial step, the anchor nodes are active, while the agent nodes 
remain silent. Then, when the inactive agent nodes are in 
connection with the active nodes and the number of these 
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active nodes satisfies the threshold value, the agent nodes 
are gradually activated by exploiting the bootstrap percolation 
strategy, which simply spreads activation from the already 
activated nodes to the hitherto inactive nodes. To elaborate 
a little further, let us define Ai as the set of agent nodes 
that are activated on layer /, where layer I is defined as the 
specific set of agents having the same level of confidence in 
position estimation. Explicitly, a coarse-grained definition of 
the confidence in the position estimate at node i is given by 


Cf{N,) 


^ 1 , Ni = l, 
< 2 , Ni = 2 , 

3, Ni>3, 


( 10 ) 


where we have Ni = It is readily observed from (10) 
that the confidence of an agent node is closely related to its 
degree of connectivity with the active nodes. Note that when 
the agent node is connected to at least three active nodes, 
its belief can be acquired as a steady unimodal distribution, 
while a bimodal distribution is obtained when the agent node is 
connected to two active nodes. If an agent node is connected to 
a single active node, the best possible estimates of its location 
constitute a ring. Therefore, using an adaptive threshold of 
connection degree is considered as the soft constraint to divide 
all sensor nodes into different layers. 


Let us assume that the agent nodes in the upper layers 
Am^ ^ = 1,...,/ — 1, have been activated subject to a 
satisfactory connection degree constraint. Then, all activated 
nodes can serve as the candidate reference nodes to facilitate 
the localization of the agent nodes on layer Ai. The set of 
candidate reference nodes given by the preceding I layers is 
defined as 


n’-= UAm,m = 0,l,...,l-l. (11) 

m 

Thus, the remaining agent nodes located above the (I — l)th 
layer are included in Bi = I\1Z\ and we have = 0. 

To assist agent node i e Bi in the localization, we select 
some neighboring nodes of the agent node i from the candidate 
reference nodes to constitute a new set as follows: 

( 12 ) 

Relying on (12), A'l is regarded as the set of the activated 
neighboring nodes of agent node i accumulated from the 
preceding I layers. Considering the nodes’ confidence in their 
position estimate, we can divide the remaining agent nodes 
having unknown positions in the set Bi into three subsets 
having different connection degree thresholds. These subsets 
are denoted by Bi^c = ^ = c}, satisfying 

Bi = UBi^c^ c = 1,2,3, where c is the connection degree 
threshold. The specific nonempty subset whose elements have 
high confidence values is preferred for constructing Ai . There¬ 
fore, Ai can be expressed as 

Ai = nrg ^mnx^Cf{Bi^c), (13) 

s.t. Bi^c ^ 

We can see from (13) that Ai is determined by the particular 
subset Bi^c which has the highest confidence in the position 



Fig. 3. Information flow among the sensor nodes on the FG of Fig. 2, 
which characterizes the scenario shown in Fig. 1(b). The circles stand for the 
variable nodes, the black squares stand for the function nodes, the dashed 
arrows denote the messages originating from the lower layers, and the solid 
arrows represent the messages emerging from the upper layers. 


estimate. In the layering process, we invoke the strategy of 
bootstrap percolation. In contrast to the conventional boot¬ 
strap percolation which uses a fixed threshold, the proposed 
bootstrap percolation strategy employs an adaptive connection 
degree threshold. 

According to the space-time hierarchical framework, agent 
nodes gradually infer their locations. Consider the example 
shown in Fig. 3, where all the agent nodes A, B, C and D 
remain silent in the initial step, while the anchor nodes 1, 2, 3, 
4 and 5 are active. The anchor nodes are collectively assigned 
to the root layer, while the agent nodes A, B, C and D are 
divided into four layers according to the connection constraints 
of ( p^ . Agent node A, which represents a class of sensor 
nodes connected with at least three reference (active) nodes, 
is activated first. Then, agent node A broadcasts its position 
to the neighboring agent nodes B and C (note that node D is 
not a neighboring node of A, as shown in Fig. 1). Afterwards, 
agent node B is activated, since it is connected to a sufficient 
number of reliable reference nodes, namely to A, 4 and 5. 
Then, agent node B estimates its own position and broadcasts 
the estimate to its neighboring node C. Subsequently, the 
connection constraints of ( p^ are rechecked, and node C can 
infer its positional distribution based on the messages sent 
from reference nodes A, B and 3. Finally, node D is activated 
by reference nodes C and 5. In contrast to the traditional meth¬ 
ods of |T^ , where the sensor nodes exchange their messages 
without control, we put emphasis on the behavior constraints 
of the sensor nodes - only the active nodes are allowed to 
estimate positions and broadcast the position information. As 
a beneficial result, an effective scheme of managing the energy 
consumption and network traffic imposed by the localization 
process is realized. According to the hierarchical approach that 
relies on the bootstrap percolation strategy described by 
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Fig. 4. The message flow on the hierarchical graph. Layer m is activated 
at time tm, and layer n is activated at time tn- The solid arrows denote the 
inter-layer messages, and the dashed arrows represent the intra-layer messages. 
For inter-layer message passing, messages always flow from upper layers to 
a lower layer. For the intra-layer message passing, such as layer n, sensor 
nodes having an insufficient number of active links can exchange messages 
for the sake of assisting each other in the position estimation, whereas sensor 
nodes having a sufficient number of active links on layer m do not exchange 
messages. 


{Tg, and ( p3] ), the agent nodes satisfying the constraint of (13) 
first estimate their positions, and then assist the other hitherto 
inactive sensor nodes to complete localization. 

In the next subsection, we will elaborate both on the 
information propagation rule and on the information fusion 
among the sensors based on the proposed hierarchical graph 
model. 


C. The NBP Positioning Algorithm Based on the Space-Time 
Hierarchical-Graph 

In this subsection, an information propagation rule based on 
the hierarchical graph model of Fig. 3 is presented first. Then, 
an NBP algorithm is developed for the space-time hierarchical- 
graph. In contrast to the standard NBP positioning algorithm 
where all sensor nodes exchange messages iteratively, we 
divide the agent nodes into multiple layers and message 
exchange is extended from the upper layers to a certain lower 
layer, as shown in Section |III-B[ Based on such a space- 
time division, it can be observed that each layer is activated 
by its upper layers during different time slots, relying on 
an initialization from the root layer and on satisfying the 
constraint in Therefore, the proposed scheme has the 
potential of realizing accurate localization during different 
time slots. More particularly, in the proposed scheme we 
have three types of messages according to the corresponding 
layers, namely the message received from the upper layers 
(MUL), the message gleaned from the same layer (MSL), 
and those emanating from the lower layers (MLL). Based 
on ([ni., the MUL has a high confidence, thus we should 
make full use of it as reliable information; the MSL serves 


Fig. 5. The belief estimation process of an agent node. The solid arrows 
denote the messages from upper layers, and the dashed arrows represent the 
messages from the same layer. Agent node i updates its belief by fusing 
the incoming messages. MUL represents messages from upper layers, MSL 
represents messages from the same layer, and MLL stands for messages in 
lower layers. Messages are exchanged among sensor nodes with the aid of 
function nodes. 


as auxiliary information, when the agent nodes are connected 
to an insufficient number of reliable reference nodes; and the 
MLL is regarded as erroneous or misleading information. In 
order to take full advantage of all valuable information, whilst 
avoiding misleading information, we define the information 
propagation rule as 

1) inter-layer message passing: messages always propagate 
from upper layers to a lower layer. 

2) intra-layer message passing: messages are exchanged 
among agent nodes, while considering the connection 
degree threshold c. 

The agent nodes connected to more than two reference 
nodes (i.e. c = 3 according to •El) are not allowed to 
exchange messages (i.e. no MSL), since they are regarded 
as misleading information. By contrast, the agent nodes con¬ 
nected to less than three reference nodes (i.e. c = 1 or 2 
according to ( p3] )) are allowed to exchange their positional 
information for the sake of assisting each other in improving 
the positioning accuracy (i.e. the MSL is required when c = 1 
or 2). 

According to the proposed information propagation rule, 
the messages fiow arranged in a controlled manner along 
the above-mentioned direction is shown in Fig. 4. Based on 
this carefully controlled information propagation, a spanning 
tree is formulated by the layering process. The anchor nodes 
constitute the root layer of the spanning tree, and the agent 
nodes satisfying the constraint of ( p^ constitute the leaf 
layers. Agent nodes at different layers are activated during 
different time slots for invoking the NBP algorithm. Provided 
that the positional information of agents on one layer becomes 
convergent after several iterations, the next layer starts to 
execute the NBP algorithm. 
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Given our information propagation rule, let us now elaborate 
on the information fusion rule invoked in NBR As shown in 
Fig. 5, the agent node i on layer I receives the MUL; when 
the agent node i is connected with an insufficient number (i.e. 
less than three) of reference nodes, it also receives the MSL. 
Then, node i updates its belief by fusing the messages received 
during time slot t as follows: 

&*(xi) oc Y[ miiixi), i G Au (14) 

u^Ti 


where the set of reference nodes is given by 


T\\jAu |-Fi|<3. 


(15) 


Note that T\ contains all the activated neighboring nodes of 
node % from the upper layers according to ^V2) . When the 
number of reference nodes is less than three, Ai contains the 
agent nodes that are in the same layer, and these nodes may 
assist node i to update its belief. Furthermore, according to 
(7), the message received from reference node u by 

node i at time slot t is expressed as 




L 




m, 


— (xu)'0L(xu,Xi)dx^, ieAuue Ti, 

" (16) 


where 6^“^(x^i) is the belief at node u in the time slot f — 1, 
representing the estimated location PDF of reference node u 
in the (t — l)th time slot. In order to apply the NBP algorithm, 
we first obtain the nonparametric expression of h\{^^ in terms 
of K weighted samples as follows p8| : 


fc = 1,2, 




(17) 


which is known as the NBP marginal. Note that x) ’ is the 
k\h sample of the belief and is the corresponding 

(k) 

weight of X- \ Similarly, the nonparametric sample-based 
expression of m^^(xi), namely the NBP message, is given 
by a Gaussian mixture message generated with the aid of K 
weighted samples, i.e. 

ueTi,k = 

(18) 

(k) 

where is the kth sample of the message arriving from 
reference node u to node i, is the corresponding weight 
of and Aui is the variance of the samples, which is 

approximated as cr^^I. 

For the sake of clarity, the detailed process of computing 
the NBP messages and beliefs (marginals) is summarized in 
Algorithm 1 fT^ . Based on these NBP messages and beliefs, 
the proposed NBP-based space-time hierarchical positioning 
procedure is summarized in Algorithm 2. 

Remark 1. Algorithm 1 provides the detailed process of com¬ 
puting the NBP messages and marginals at time slot t, where 
> 1 is a sampling parameter. The parameter h is adjustable. 
A lager h indicates relying on more fused information and on 
increased processing time. 


Remark 2. By exploiting a heuristic search relying on the 
bootstrap percolation strategy, the agent nodes are divided 


Algorithm 1 Computation of NBP Messages and Marginals 
I: Compute NBP messages: Given K weighted samples 
from construct an approximation to 




(xi) from each activated neighbor s G <5^: 


2: if pcz(xs,Xi) = 1 then approximate mT(xi) with a Gaussian 
mixture: 

3: Draw random values for ^ f/[0,27r) and ~ pn 

where pn denotes the PDF of the ranging noise. 

4: Means: = xi^^ -h {dsi + n^^^)[sin(^^^^); cos(^^^^)] 

5: 


Weights: = 


(fc) 




6: Variance: As* = K~^ • Cov[m^^^] 

7: end if 

8: Compute NBP marginals: Given several Gaussian mixture 
messages , Aui}, draw samples from 

^n“^(Xn): 


9: for each observed neighbor u G do 
10: Draw samples {x-^^} from each message 

11: Weight by wf ^ V EnsjC; 

12: end for 

13: From these hK locations, resample M times according to 
weights (with replacement) to produce K equal-weight samples. 


Algorithm 2 The Space-Time Hierarchical Positioning Algo¬ 
rithm_ 

I: Initialization: The root layer is = J, and the reference nodes 
set is 72° = J. 

2: for / = 1 to L do {iteration index} 

3: Determine sensors to be located on layer / Select 

the agent nodes in Ai using GD- 
4: for i ^ Ai do (the reference nodes set of node i is the set 

= =Sin'R.‘} 

5: if Cf(Ai) = 3 then 

6: Intra-layer messages are not in use 

7: else if Cf{Ai) < 2 then 

8: Add assistant reference nodes from the intra-layer 

nodes using ^ U Ai 

9: end if 

10: end for 

II: end for 

12: Information fusion: Sensors in Ai receive messages from the 
preceding I layers TZ\ and update their beliefs. 

13: for z G do (start the NBP computation using Algorithm 1 in 
parallel} 

14: for t = 1 : T do 

15: Compute NBP messages from its reference nodes Ti 

16: Compute NBP marginals for agent i 

17: end for 

18: end for 

19: Reference nodes updating: Update the reference nodes set 
using ^ yj Ai 


into multiple layers. Thus, a space division of the positions is 
obtained, which is capable of mitigating the infiuence of loops 
in a graph. Based on the proposed information propagation 
rule, the messages fiow in controlled directions and the MLL is 
regarded as misleading information, which is censored. Addi¬ 
tionally, localization estimation is implemented in a layer-by- 
layer manner, where the agents gradually infer their positions 
by selecting reference nodes from the upper layers. In this 
way, the network traffic overhead can be dramatically reduced. 
Moreover, for the agent nodes that have an insufficient number 
of reliable reference nodes, the messages can be exchanged 
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TABLE II 

Complexity Comparison in Terms oe the Average Number oe 
Links Used in the Entire Network 


Algorithm 

NBP 

NBP-BFS Hierarchical NBP 

Complexity 

7rV2p2 

2N -2 


among agents on the same layer (intra-layer message passing). 
As a result, all helpful information can be exploited for 
improving the localization accuracy. 

IV. Computational Complexity and Network 
Traffic Overhead 

In this section, we analyze both the computational com¬ 
plexity and network traffic overhead of the proposed NBP- 
based space-time hierarchical positioning algorithm. Both the 
computational complexity and the network traffic are related 
to the number of actively participating links on the graph, 
hence we adopt the average number of links used in the entire 
network as an indirect metric for comparing the complexity of 
different NBP algorithms. In particular, a comparison is made 
among the proposed algorithm, the standard NBP algorithm 
and the low-complexity NBP defined over spanning trees with 
breadth first search (NBP-BFS) pT| . We consider M anchor 
nodes placed at fixed positions and N agent nodes randomly 
distributed over an a x a square area. The communication range 
is given by R = pa^ where p is a proportionality coefficient. 
Then TVmsg. which is the average number of messages fused 
by an agent node, can be written as 

Nmsg = ^ ^ 7r(pa)^ = Nnp^. (19) 

In the standard NBP algorithm, all messages are employed for 
information fusion. Thus, the computational complexity of the 
standard NBP algorithm is characterized by 

OnBP = N • Nrnsg = (20) 

When using the NBP-BFS algorithm, — 1 edges are retained 
in the graph and there are two messages on each edge, thus 
the computational complexity of the NBP-BFS is determined 
by 

Obfs = 2{N-1). (21) 


By comparison, as far as the proposed hierarchical NBP 
scheme is concerned, all the sensor nodes are divided into L 
layers. Therefore, the computational complexity of the pro¬ 
posed algorithm in each iteration is approximately determined 
by 


^layer — A/^ ’ Allayer — ^ 


N A^V 


L 




Onbp 


( 22 ) 

where N stands for the average number of agent nodes in 
each layer, Mayer represents the average number of messages 
that have to be fused. In particular, if L = 1, the proposed 
algorithm shares the same complexity as the standard NBP. 
The complexity comparison results of the three algorithms are 
summarized in Table [III 
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Fig. 6. Complexity comparison of different NBP algorithms in terms of the 
average number of links used in the whole network, where p = 0.2 is fixed 
and N changes from 50 to 100. 



p [%] 

Fig. 7. Complexity comparison of different NBP algorithms in terms of the 
average number of links used in the whole network, where W = 100 is fixed 
and p changes from 0.2 to 0.6. 


If we have Oiayer < Obfs, the following conditions should 
be satisfied, 

^< 2 (iV-l), ^ < 2, and JVV < 

which implies that the proposed algorithm tends to exhibit a 
reduced complexity when more layers are created. As shown 
in Fig. 6, we consider a typical large-scale network (H having 
13 anchors distributed in a square area of 100m x 100m, i.e. 
a = 100, and the transmission radius is fixed io R = 20m, 
hence p = 0.2 and N changes from 50 to 100. It can 
be observed from Fig. 6 that when the number of agent 
nodes increases, the complexity of the standard NBP also 
increases, since there are more messages that have to be 
fused, while the proposed algorithm dramatically reduces the 
complexity. Fig. 7 shows our complexity comparison of the 
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Fig. 8. Deployment of anchor nodes. The circles indicate the radio range of 
the anchors having a transmission radius of = 12m. 

different NBP algorithms subject to the fixed N = 100 and 
to the transmission radius varying from 20m to 60m, i.e. p 
ranges from 0.2 to 0.6. We can see that upon increasing 
the transmission radius, the average number of links in the 
network increases, hence resulting in an increased complexity 
for the NBP. By comparison, the proposed hierarchical scheme 
significantly reduces the number of participating links, thus it 
has a reduced complexity. Additionally, we can see from Fig. 7 
that when there are less layers, the proposed algorithm suffers 
from a higher complexity than NBP-BFS. Moreover, when the 
number of layers L increases, the complexity of the proposed 
algorithm decreases and it has the potential to become lower 
than that of NBP-BFS. 

V. Simulation Results and Discussions 

In this section, numerical simulations are carried out for 
evaluating the performance of the different positioning algo¬ 
rithms considered. In summary, compared to the standard NBP 
algorithm, the proposed space-time hierarchical localization 
scheme attains a performance gain by censoring the MLL 
that contains misleading information. In contrast to the low- 
complexity NBP-BFS of pT| , which runs the NBP on a 
spanning tree without considering the appropriate number of 
reference nodes, the proposed scheme runs NBP on a space- 
time hierarchical-graph, which ensures that most of the agent 
nodes share a sufficiently high number of reference nodes for 
localization, thus a performance gain is achieved. 

A. Simulation Setup 

We consider a 50mx 50m square area in our simulations, 
where a number of anchor nodes are placed at fixed positions 


TABLE III 

Deployments oe dieeerent networks 



Network 1 

Network 2 

Network 3 

Agents 

13 

13 

9 

Anchors 

100 

50 

100 


and a large number of agent nodes are randomly distributed 
across the area. The network deployment parameters are 
summarized in Table III, and the a priori positions of the 
anchors are shown in Fig. where denotes the loca¬ 
tions of the anchor nodes in Network 1 and Network 2; 
denotes the locations of the anchors in Network 3; and 
denotes the locations of the agents. The range measurement 
noise is modeled using the Gaussian distribution with zero 
mean and a standard deviation linearly depending on the 
distance. Thus the estimated distance is given by d = d + 

Us ~ A/’(0, (cTo + /Co- X d)^), where cfq and ka- are constant 
values. In the simulations, it is assumed that ao = 0.2m 
and ka- = 0.01. The transmission radius is 12m, and we use 
K = 200 samples for approximating the beliefs and messages. 
The central processing unit (CPU) utilized for the simulations 
is an Intel Core 15-4570 CPU. Two criteria are utilized for 
evaluating the performance, including the CPU running tim^ 
and the error distribution CDF, which is given by 

CDF(e,^) =E{I{||x-x|| < e,4}, (23) 

where x = [xi,..., xat], x = [xi,..., xat], and eth is the 
positioning error threshold of an agent. In the simulations, the 
estimated parameters are exchanged among the sensor nodes 
10 times. All the positioning errors are obtained for 1000 
randomly generated networks. 

B. Performance Results and Discussions 

Fig. 9, Fig. 10 and Fig. 11 show the positioning performance 
of the proposed NBP based space-time hierarchical positioning 
algorithm subject to differently initialized connection degree 
thresholds given in •ini' for the three networks specified 
in Table III, respectively. The positioning performance of 
Network 1 is better than that of the other two networks, 
because Network 1 has enough anchor nodes and dense agents. 
Note that when the connection degree threshold is initialized 
to zero, the proposed space-time hierarchical NBP reduces to 
the standard NBP. It can be seen that the proposed space- 
time hierarchical NBP algorithm using a connection degree 
threshold higher than two outperforms the standard NBP 
algorithm (i.e. the connection degree threshold is initialized 
as c = 0), due to its resilience against error propagation. 
Additionally, we evaluated the scenario where different hi¬ 
erarchical graphs are generated corresponding to differently 
initialized connection degree thresholds. We observe from 
Fig. 9, Fig. 10 and Fig. 11 that there are no significant 

^ Apart from the metric of “the average number of links used in the entire 
network”, which is used in Table II, the CPU running time is invoked in 
the simulations as a complementary metric for numerically characterizing the 
computational complexity. 
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Fig. 9. Comparison of positioning error CDFs of hierarchical NBP under 
differently initialized connection thresholds for Network 1 specified in Table 
III. 



Fig. 10. Comparison of positioning error CDFs of hierarchical NBP under 
differently initialized connection thresholds for Network 2 specified in Table 

III. 

positioning performance differences for the proposed space- 
time hierarchical NBP algorithm, when the connection degree 
thresholds are increased from zero to two. Moreover, the 
space-time hierarchical NBP algorithm exhibits a significantly 
improved performance when the connection degree threshold 
is higher than two. Fig. 12 shows that the CPU running time 
of the proposed space-time hierarchical NBP algorithm can 
be reduced when more layers are created, which is consistent 
with the complexity analysis results shown in Section IV. 

The CDF comparison of different NBP algorithms for 
Network 1 and Network 2 is shown in Fig. 13, where it 
is observed that the CDF of the NBP-BFS | [3T| indicates a 
poor positioning performance, because insufficient information 
is gleaned. By comparison, the space-time hierarchical NBP 
(where the connection degree threshold is initialized as three) 



Fig. 11. Comparison of positioning error CDFs of hierarchical NBP under 
differently initialized connection thresholds for Network 3 of Table III. 



Fig. 12. CPU running time comparison of hierarchical NBP under different 
connection thresholds for the three networks of Table III, respectively. 

makes full use of all beneficial information and mitigates 
error propagation, thus attaining an improved positioning 
error performance. We can see that the proposed space-time 
hierarchical NBP shows the best performance, and most of the 
sensor nodes are located within a positioning error of less than 
Im. In addition to its beneficial performance gain, we can see 
from Table IV that the proposed space-time hierarchical NBP 
algorithm significantly reduces the complexity, hence incurring 
a reduced CPU running time. We observe that with the aid of 
the proposed space-time hierarchical strategy, the NBP can be 
executed roughly 15 and 10 times faster than the standard NBP 
for Network 1 and Network 2, respectively. 

VI. Conclusions 

In this paper, a space-time hierarchical-graph was proposed 
for efficient cooperative positioning in WSNs. Specifically, 
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TABLE IV 

CPU RUNNING TIME COMPARISON OF DIFFERENT POSITIONING ALGORITHMS FOR NETWORK 1 AND NETWORK 2. 


Positioning algorithms 

NBP |18), |29 

NBP-BFS 

NBP-MIN jE 

Hierarchical NBP 

CPU running time of Network 1 [s] 

71.03 

10.53 

15.24 

4.80 

CPU running time of Network 2 [s] 

15.01 

4.59 

5.94 

1.47 



Fig. 13. CDF comparison of different NBP algorithms for Network 1 and 
Network 2 of Table III. 

a bootstrap percolation strategy relying on soft constraints 
was utilized for controlling the actions of sensor nodes. As 
a benefit of the hierarchical mechanism, agent nodes that 
satisfy our connection constraints are activated gradually for 
estimating their positions, while the sensors that do not 
satisfy the connection constraints remain silent. Hence, the 
energy consumption/network traffic overhead imposed by the 
localization process is significantly reduced compared to that 
of the benchmarking schemes considered. Additionally, an 
information propagation rule was designed by exploiting the 
hierarchical graph. As a result, the messages always fiow in the 
direction from the sensor nodes of higher confidence to those 
having a lower confidence. Therefore, the error propagation 
can be effectively mitigated, yielding an accelerated position¬ 
ing convergence. In our future work, the bootstrap percolation 
strategy relying on soft constraints and other features of 
WSNs (e.g. multi-hop) may be exploited in network forming, 
which is expected to provide an improved performance gain 
or to reduce the network traffic overhead imposed by the 
localization process. 
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